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ABSTRACT 

We present long-slit, optical spectra of the merging system Arp 220, taken with the 
William Herschel Telescope (WHT) on La Palma. These data were taken as a part of a 
large study of ultraluminous infrared galaxies (ULIRGs) with the aim of investigating 
the evolution and star formation histories of such objects. Spectral synthesis modelling 
has been used to estimate the ages of the stellar populations found in the diffuse light 
sampled by the spectra. As the closest ULIRG, it proved possible to perform a detailed 
study of the stellar populations over the entire body of the object. The data show 
a remarkable uniformity in the stellar populations across the full 65 arcsec covered 
by our slit positions, sampling the measurable extent of the galaxy. The results are 
consistent with a dominant intermediate-age stellar population (ISP) with age 0.5 < 
t/sp < 0.9 Gyr that is present at all locations, with varying contributions from a 
young (< 0.1 Gyr) stellar population (YSP) component. However, it is notable that 
while the flux contribution of the YSP component in the extended regions is relatively 
small (< 40%), adequate fits in the nuclear region are only found for combinations 
with a significant contribution of a YSP component (22 - 63% ). Moreover, while 
a low intrinsic reddening {E{B — V)< 0.3) is found for the ISPs in the extended 
regions, intrinsic reddening values as high as E{B — V) ^ 1.0 are required in the 
galactic center. This clearly reflects the presence of a reddening gradient, with higher 
concentrations of gas and dust towards the nuclear regions, coinciding with dust lanes 
in the HST images. Overall, our results are consistent with models that predict an 
epoch of enhanced star formation coinciding with the first pass of the merging nuclei 
(represented by the ISP), with a further episode of star formation occurring as the 
nuclei finally merge together (represented by the YSP and ULIRG). 
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1 INTRODUCTION 

Since the early 1980s, several studies have revealed the 
presence of galaxies with the spectral ene rgy di s tribu- 



tion domina te d by i nfrared emissio n ( Houck fc et al j_ 1984 



Houck et al.l. Il985l: ISoifer fc et all . Il984a .h: ISoifer et al 



19871 : iLe Floc'h et al.l . 12005^ . These galaxies are classified 
Lq) or Ultraluminous [Lir > 10^^ 
Their prodigious 



as Luminous (Li,- > 10 
L0) infrared galaxies (LlRGs/ULIRGs 
infrared emission is generally attributed to the optical/UV 
light of luminous central sources reprocessed by dust. Star- 
burst activity is ongoing in most, if not all, of these sources , 
coexisting with AGN activity in some cases jGenzel et al.l . 
1 19981 : ISurace fc Sandersl Il999l : IVeilleux et all . 1 19991 ). There- 
fore these objects provide us with an excellent oportunity 
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to study both AGN and starburst phenomena. Moreover 
ULIRGs are almost unv ariably associated wi t h gala xy merg- 
ers and interactions fsee lSanders fc Mirabell l|l996l) for a re- 
vision). Thus, they also represent ideal objects to test models 
of galaxy evolution via major galaxy mergers, (e.g. Mihos fc 
Hernquist, 1996; Barnes & Hernquist, 1996; Springel et al., 
2005.) 

In this context Arp220, as the closest ULIRG (z — 
0.018), is a key object for our understanding of star- 
formation in galaxy mergers. At inf rared waveleng t hs, it 
appears as a double nucleus system (jGraham et al.l . Il990l : 
IScoville et al] . 119981 ) separated 0.95 arcsec, with two tidal 
tails visible at optical wavelengths ( Josep h fc Wrightl . ll985l ). 
Although the mid-IR spectra suggest starburst activity 
as the main s ource of power for the infrared continuum 
iGenzel et al.l . 1 19981 : iLutz et all , 
sified as a LINER in the optical 



19991). the galaxy is clas- 



Veilleux et al] . 1 19991 ). Arp 
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Figure 1. Hubble space Telescope (HST) images taken with the F814W filter, showing the locations of the slits and the extraction 
apertures. Left: WFPC2 image showing both the nuclear region and the extended regions. Right: ACS HRC Image showing the 
region within the dashed-line box in the left pannel. 



220 was also observed by IScoville et all at 1.3 and 

2.6 /xm, sampling the CO (2-1), CO (1-0) and dust contin- 
uum emission. They determine a total H2 mass of 3.2 x 
10^° M0, of which 2/3 is contained within 250 pc and con- 
fined to a thin, centrally located, dense disk, embedded in 
a CO disk extended by ~1 Kpc . The dynamics of the nu- 
clear region are consistent with a collision between prograde 
and retrograde disks overlapping with the main gas disk ( 
Mundel et al . 2001 , see figure 7 in their paper). Recently, 
IWilson et al.l l|2006l l performed a detailed photometric study 
of the star clusters in Arp220. They found evidence that the 
cluster population can be divided into two groups: a cen- 
trally located young population of age ty sp < < 10 Myr, 
and an intermediate-age population of age ~300 Myr spread 
towards the north of the galaxy. However, given the photo- 
metric uncertainties, it is not straightforward to break the 
age/reddening degeneracy problem unambiguously in such 
studies, based only on colour-colour diagrams with relatively 
few filters. 

To date, there have been few spectroscopic studies of 
stellar populations in ULIRGs. To remedy this situation we 
present in this paper a detailed study of the stellar popula- 
tions in Arp220, based on spectroscopic observations of the 
extended diffuse light along three slit positions. The results 
are discussed in the context of evolutionary models of star 
formation in mergers. 

Throughout the paper we assume Ho = 71 km s~^, 
f2o = 0.27 and Q.a = 0.73, resulting in a scale of 0.363 kpc 
arcsec"^ and a distance of 77.6 Mpc at z = 0.018. 



2 OBSERVATIONS AND DATA REDUCTION 

Spectra of Arp220 were taken in July 2005 and 2006 with 
the ISIS dual-beam spectrograph on the 4.2-m William Her- 
schel Telescope (WHT) on La Palma. The R300B grating 



with the EEV12 CCD, and R316R grating with MARC0NI2 
CCD were used on the blue and the red arms respectively. 
A dichroic cutting at 5300A was used during the observa- 
tions, leading to a useful wavelength range in the rest frame 
of 3000 - 5OOOA in the blue, and 5000 - 7800A in the red. 
The spectra were taken along three slits positions: PA 160, 
PA 75 and PA 75* (10 arsec offset to the north of PA 75), 
covering the nuclear region and also the tidal tail observed 
towards the North- West of the galaxy. Details of the obser- 
vations are given in Table 1. Figure 1 shows the slit posi- 
tions superimposed on HST images taken with the WFPC2 
(HST proposal 6346, PI K. Borne) and the ACS HRC (HST 
proposal 9396, PI C.Wilson) camera. To find the precise lo- 
cations of the slits on the images we first convolved the im- 
ages with Gaussian profiles to simulate the seeing conditions 
during our WHT observations. Several spatial profiles were 
then extracted from the images and compared with those 
of spatial slices extracted from the spectra using a wave- 
length range as close as posible to that of the images, until 
a match was found. All exposures were taken with the slit 
aligned along the parallactic angle, in order to minimize the 
effects of the differential refraction. 

The data were reduced (bias subtraction, flat field cor- 
rection, cosmic ray removal, wavelength calibration and fiux 
calibration) and straightened before the extraction of the in- 
dividual spectra using the standard packages in IRAF and 
the STARLINK packages FIGARO and DIPSO. The wave- 
length calibration accuracy, calculated a s the mean value 
of th e difference between the published (jOsterbrock et aP . 
Il996l l and the actual wavelengths of the night sky lines was 
estimated as < 0.34A and < 0.45A in the blue and the 
red respectively. The spectral resolution, calculated as the 
mean value of the sky line width (FWHM) was 5.8OA m 
the blue and 5.25A in the red. The relative flux calibration 
accuracy, based on the observation of several spectrophoto- 
metric standard stars, is estimated as ± 5% over the entire 
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Figure 2. The Balmer break ratio Rs plotted as a function of 
position along PA 160. The small variations of the ratio suggest 
a remarkable uniformity for the stellar populations along the ex- 
tension covered by the slit position. The centroid of Ap D is used 
as the reference point. Positive x = North, Negative x = South. 



spectral range. This is confirmed by the excellent match be- 
tween the blue and the red spectra for all the apertures ex- 
tracted. Prior to the modelling, t he spectra wer e corrected 
for Galactic extinction using the ISeatonl ll 19791) reddening 
law an d the value of E{B ~V) = 0.051 from lSchlegel et all 
(| 19981 ). 



3 RESULTS 

In order to gain a preliminary impression of the variations 
of the stellar populations, in Figure 2 we plot the flux ratio: 



Rl 



(■4080 

J4OOO ^^'^^ 
(■3700 



(1) 



as a function of position along PA 160. This ratio, which is 
designed to provide a sensitive measure of the Balmer break, 
shows remarkably little variation across the galaxy (typical 
variation ± 10%) except in the region of the nucleus (within 
± 5" of the center) which is subject to enhanced reddening. 

To further investigate the stellar populations, 20 aper- 
tures were extracted along the three slits positions with the 
aim of adequately sampling the different regions observed in 
both the 2-D spectra and the available HST images. These 
apertures are labelled in Figure 1. Three larger apertures, 
including all the small apertures for each slit PA, were also 
extracted, the latter are labelled as KVtotal in Fig 3 and 
Table 2 and 3. The apertures extracted for the three slit po- 
sitions are shown superimosed on spatial profiles of the 2-D 
spectra for the wavelength range 4400 - 4600A in Figure 
3. The apertures were chosen to be wide enough to have a 
sufficiently high S/N ratio for further analysis: S/N > 15. 
We would like to emphasize that the S/N is considerably 
higher than this for most of the apertures extracted from 
the different regions of the galaxy. A sample of extracted 
spectra is shown in Figure 4. The strong Balmer lines and 
Balmer break present across the full measureable extent of 
the galaxy (65 arcsec) indicate that young and/or interme- 





-40 -30 -20 



Arcsec 



Figure 3. Spatial profiles of the 2-D spectra for the wavelength 
range 4400 — 4600A showing the positions of the extraction aper- 
tures. The centroid of Ap D is chosen as the reference point for 
the three slit position, PA160, PA75 and PA75*. PA 160: Posi- 
tive X = North, Negative x = South. PA75 and 75*: Positive x = 
West, Negative x = East. 



diate age stellar populations are present in all spatial aper- 
tures. The overall SEDs also show signs of enhanced redden- 
ing in the nuclear regions (apertures D and E, PA 160). 

In order to obtain a more detailed information about the 
stellar populations we have modelled the spectra extracted 
fro m the apertures using th e evolutionary synt hesis models 
of iGonzalez Deleado et aP (|2005,) and .Bruzual fc CharlotI 
(|2003h . These models incorporate high and intermediate 
spectral resolution stellar libraries. Such resolutions are re- 
quired to isolate the nebular and the stellar absorption con- 
tributions in the spectra of this type of galaxies. To generate 
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Date 


Arm 


set up 


Slit PA 


Slit width 


Exposure Time 


Seeing 






(CCD + gratting) 




(arcsec) 


(sec) 


(arcsec) 


July 2005 


B 


EEV12 + R300B 


160 


1.5 


900 


1.85 


July 2005 


R 


MARCONI2 + R316R 


160 


1.5 


900 


1.85 


July 2006 


B 


EEV12 + R300B 


75 


1.5 


900 


0.65 


July 2006 


R 


MARCONI2 + R316R 


75 


1.5 


900 


0.65 


July 2006 


B 


EEV12 + R300B 


75* 


1.5 


1200 


0.75 


July 2006 


R 


MARCONI2 + R316R 


75* 


1.5 


1200 


0.75 



Table 1. Summary of the spectroscopic observations of Arp220. 



these templates we have assumed a ISalpeted (|l955l ) initial 
mass functio n, instantan e ous s tarburst and solar metallic- 
ity (sec Tad hunter et al.l . |2005| . for a justification of these 
assumpti on). We then reddene d the synthetic spectra using 
both the ICalzetti et al.l (120001) reddening l aw, appropriate 

ie|s ■ 



for starburst galaxies, and the lSeatonI (Il979l ) reddening law 
representing the Galactic extinction case. Overall, there is 
good consistency between results obtained with the two red- 
dening laws, demonstrating that the main results are not 
sensitive to the details of the reddening law assumed, at 
optical wavelengths. 

To determine the detailed properties of the young stellar 
populations we have to fit the whole observed spectrum. We 
have used two different approaches to perform the fit and de- 
termine the best models results. They involve STAR LIGHlH 
Cid Fcrnandcs ot al, 2005) and CONFIT ( Robinsoji_eLaI 
iTadhunter et all . 120051 : iRodn'guez Zaurm et al 



2007 



Holt et al.l . I2OO7I ) codes. The former allows up to 11 stellar 



components but with the same extinction for all compo- 
nents. On the other hand, the CONFIT approach assumes 
only two stellar population components of different age, but 
allows these components to have different reddening. More 
details of the two techniques are given in sections 3.1 and 
3.2. 

Throughout this paper, we define young stellar popu- 
lations (YSPs) as stellar components with ages tysp < 0.1 
Gyr, intermediate-age stellar populations (ISPs) as stellar 
components with ages in the range of 0.1 < tjsp < 2 Gyr and 
old stellar populations (OSP) as a component with age 12.5 
Gyr. Note that we have not used stellar populations with 
ages in the range of 2 Gyr < tsp < 12.5 Gyr. Generally, for 
models that include young and intermediate/old stellar pop- 
ulations, which is the case for both CONFIT (2 stellar com- 
ponents plus a power law in some cases) and STARLIGHT 
(11 stellar components), it is not possible to distinguish be- 
tween ages within this range, and therefore we decided not 
to use such stellar populations during the modelling analysis 
described below. 



3.1 Results from STARLIGHT 

STARLIGHT fits simultaneously the continuum points and 
the absorption lines (high order Balmer lines. He lines. 
Call H, Call K, G-band...). The program reads an input 
file that labels the spectral windows that contain emission 
lines, and these spectral ranges are excluded of the fitting. 



www.starlight.ufsc.br 



In these models, we have only used the SSP from the evo- 
lutionary modelfl The intrinsic extinction is modelled by a 
foreground dust screen, and parametrized by the V-band ex- 
tinction, Av- The output is the population vector that rep- 
resents the fractional contribution of each SSP of a given age 
and metallicity to the model flux at a given wavelength. The 
fit is carrie d out with a simulate d annealing plus Metropo- 
lis scheme (|Cid Fernandes et al.l . 12005 ) , which searches for 
the minimum chi-squared between the observations and the 
combined models. 

For this work we used 11 SSP, corresponding to ages 4, 
5, 10, 25, 40, 100, 280, 500, 890 Myr, and 1.25 and 14 Gyr 
from models bv lGonzalez Delgado et al] (|2005| ). We use so- 
lar metallicity, and the Calzetti et al. (2000) extinction law. 
The results obtained are shown in Table 2 and are summa- 
rized as follows: 

• OSP contribute very little to the optical continuum. 
Only a few percent of the flux at 4020A is due to stellar 
populations older than 2 Gyr. 

• ISPs with ages in the range of 0.5 Gyr < tjgp < 0.9 
Gyr account for > 70% of the optical continuum in all aper- 
tures apart from apertures D and E in PA 160, sampling the 
nucleus of the galaxy. 

• The contribution of YSP (< 100 Myr) is signiflcant (> 
30%) in the nuclear regions. 

• While low intrinsic reddening is found in the extended 
regions (0.0 < E{B ~ V) < 0.4 ), large intrinsic reddening, 
E{B — V) ~ 0.7 - 0.8 is required towards the centre of the 
galaxy. 

• The modelling results found for the three large aper- 
tures are consistent with those of the smaller ones and with 
each other, as shown in the table. 

Similar conclusions are obtained if we use 
iBruzual fc CharlotI l|2003l ) synthetic SEDs. 

A drawback of the STARLIGHT code is that all steUar 
population components have the same reddening. Therefore, 
we have also used the CONFIT code, allowing the two stellar 
components used in this approach to have different redden- 
ings. 

3.2 Results from CONFIT 

The CONFIT approach consists of a direct flt of the 
overall continuum shape of the extracted spectra us- 

^ Starlight can also add a power law to the SSP to fit the con- 
tinuum, but these models do not consider any contribution by a 
power law 
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ing on a minimum techn ique l|Tadhunter et all l2005l : 
iRodn'guez Zaun'n et all . 120071 ') . CONFIT is based on a "sim- 
plest model" approach, i.e. we fit the minimum number of 
stellar components required to adequately model the data. 
The main advantage of this approach is that it allows a wide 
range of reddening values for each of the components to be 
explored. For this work we have used different combinations 
of two stellar components (ISP +YSP) with a range of red- 
dening for each component. For each spectroscopic aperture 
the continuum flux was measured in several wavelength bins 
(~70) chosen to be as evenly distributed in wavelength as 
possible, and to avoid strong emission lines and atmospheric 
absorption feature^ The model flt results are quantifled in 
terms of the percentage contribution of the different stellar 
components in a normalising bin deflned as 4700 - 4800 A. 
Due to the uncertainties associated with the flux calibra- 
tion, as well as the synthetic templates themselves, it is not 
correct to constrain the total flux in the normalising bin to 
be exactly the integrated flux of the observed spectrum in 
the bin. Therefore the CONFIT codes allows the model flux 
in the normalising bin to vary up to 125% of the measured 
flux in that bin. 

Potentially adequate fits to the SEDs (reduced < 1.0) 
can be found for combinations of small contributions of 
highly reddened YSP plus a very old stellar population, or 
large contributions of low reddened ISP with an old stel- 
lar component contributing much less to the overall stellar 
population in the aperture (i.e. the age/reddening degener- 
acy problem). To break this degeneracy we have examined 
detailed fits to the spectra, considering as valid those mod- 
els that fit the important absorption features of the spec- 
tra adeq uately (e.g. Call K, G-band, h i gher order B aime r 
lines: seelConzalez Delgado et all (|200ll ): lHolt et al.l l|2003ll : 
iRodrfguez Zaun'n et al.1 l|2007h '). 

Table 3 summarizes the results obtained with the CON- 
FIT code. Adequate fits for all 23 apertures extracted from 
the three slit positions are found with models combining 
two stellar components, an ISP and a YSP, with no need of 
an OSP. Figure 4 shows example fits for apertures selected 
with the aim of sampling a range of regions of the galaxy. 
Note that apertures Ap D and E for PA 160 sample the cen- 
tral 2.5 kpc diameter region, corresponding to dust lanes in 
the HST images. It is clear from the figure that the spectra 
extracted for these apertures are highly reddened. As well 
as the best-fitting models, detailed fits in the region of the 
Balmer lines (3700 - 4300A) are also shown in Figure 4. 

Overall, the results reveal a remarkable uniformity in 
the estimated stellar ages, with an ISP (0.5 Gyr < tisp < 0.9 
Gyr) dominanting the optical emission in all apertures, and 
a YSP (< 0.1 Gyr) making a relatively small contribution (< 
40%) in most apertures. However note that in the central 2.5 
kpc of the galaxy (Ap D and E, PA 160), adequate fits are 
only found for combinations with a significant contribution 
of a YSP component (22 - 63% ). In terms of reddening, 
while the results found in the extended regions are consistent 
with low intrinsic reddening, one or both stellar components 



^ Note that it is a further difference with the STARLIGHT code 
which fits every pixel apart from regions (e.g. those containing 
emission fines or atmospheric emission features) that are specifi- 
cafiy masked out. 



must be highly reddened in the central 5 kpc of the galaxy, 
clearly indicating high concentrations of gas and dust in the 
nuclear regions. 

At this stage, it is important to add a caveat about the 
spatial extent of the YSP component. In some of the ex- 
tended regions of the galaxy (Ap E and F, PA 75*; Ap A, F 
and G PA 75; Ap A, PA 160), good fits are also found with a 
minimal contribution of a YSP template. This is consistent 
with the STARLIGHT resuffs, where the ISP contribution 
found for these apertures is > 90 % of the total fiux. On the 
other hand, a YSP component is always required to model 
the data for the apertures sampling the center of the galaxy. 
Therefore, while the presence of a YSP in the extended re- 
gions remains uncertain, the presence of such a component 
in the galactic center is ess ential. This re s ult is consistent 
with the age of 10-100 Myr (iMundefl et al.l . 120011 1 estimated 
for the duratio n of the starbur s t asso ciated with the galac- 
tic superwind l|Heckman et al.l . Il996l ) in the central region 
of the galaxy, and also wi t h the merger simulations predic- 
tions llMihos fc HernqulsR 1 19961 : [Barnes fc Hernguistl [l996l : 



ISpringel et al.l , l2005l ) for starbursts characteristic of the final 
stages of major gas-rich mergers. 

Again, the results for the three large apertures show 
remarkable consistency, both with each other, and with the 
smaller aperture results presented in the tables. 

Overall, the results obtained with CONFIT are broadly 
consistent with those of STARLIGHT approach, and clearly 
reinforce the idea that an ISP dominates the optical emission 
in Arp220. 



3.3 The Hidden Starburst 

Our analysis has so far concentrated on the optically visible 
star formation. It is interesting to consider the extent to 
which this represents the full extent of star formation in 
Arp220. 

Previous studies based on radio observations 
(Smith et al.l. Il998f) or the galactic superwind 



I bmith et al.l. litJijal) or tne galactic superwind 
( Heckman et al.l . Il996l : iMundeU et"all 120011 ) . have found 
star formation rates (SFR) up to 100 M© yr"^ for Arp 220. 
An independent estimate of the SFR can be obtained using 
the publis hed ffux of the Bra l ine, F(Bra) = 1.8 x 10"" erg 



(ICenzel et al. J 1998) . Using t he IR extincti o n laws 

of llndebetouw et all (|2005i ) and Rieke fc Lebofskvl (ll985D 
along with the estimated visual extin ction for the nuclear 
regi on of Arp 220 (A„ ~ 45 mag, ISturm et all (|l996l ): 
Gen zel et al.l (|l998l )). we derive a Bra luminosity of L(Bra) 



theory we then derive an Ha luminosity of L(Ha) = 3.3 x 
10^^ erg s-^ It is now straightforward to estimate the SFR 
in the galaxy using the relation ship between Ha luminosity 
and SFR given by iKennicutti (il99§ ). We find a SFR ~ 
260 Mq yr~^, enough to power the entire mid- to far-IR 
luminosity observed in Arp220. Moreover, the absence of 
mid-IR, high-ionization emission lines ([NeV] A14.322//m, 
OIV]A25.890/im. Genzel et al., 1998) suggests no hidden 
AGN actvity, reinforcing the idea that starburst activity 
powers the ULIRG Arp220. In such a case, it is also possible 
to measure the S FR of the galaxy using the IR-luminosity 
(iKennicuttll 19981). For a val ue of Lir = 1.3 x 10^^ L© for 
Arp 220 JSoifer et al.l , Il987l ) we obtain a SFR = 224 Mq 
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yr~^ which is remarkably consistent with the value found 
using the Bra flux. 

However, there is no clear evidence for the pres- 
ence of ongoing star formation activity at optical wave- 
lengths, based o n HII region-like optical emission line ratios 
l|Veilleux et al.l , Il999l ). In order to measure the bolometric 
luminosity associated with the visible stellar populations de- 
tected in the optical, we added the spectra extracted for all 
the apertures of the region limited by the dashed-line box 
in Figure 1. Since our spectroscopic slits used do not cover 
the entire extent of the central region, a large aperture cho- 
sen to match the extent of the box was extracted from an 
HST ACS (F435W) image (HST proposal 9396, PI C. Wil- 
son). We then scaled the combined spectrum to match the 
flux measured from the aperture, in order to account for 
possible flux losses. Assumming an ISP -I- YSP combination 
consistent with the modelling results of section 3. and cor- 
recting for reddening effects, we determined the bolometric 
luminosity of all the components in the scaled spectrum and 
added them, obtaining a total bolometric luminosity Ltoi '~ 
1.6 X 1O"L0. This is an order of magnitude less than the 
mid- to far-IR luminosity of the source. Therefore, especially 
considering that not all the optical emission is absorbed by 
dust, it is clear that the stellar populations detected at op- 
tical wavelengths cannot power the far-IR luminosity of the 
source; most of the ongoing star formation in the nuclear 
regions is hidden by dust. 

3.4 Summary of the key results 

To summarize tha main results, the modelling of our spectra 
provide evidence for three stellar populations in Arp220: 

• Dominant ISP. We have found an intermediate-age 
stellar population, with ages 0.3 Gyr < t/sp < 0.9 Gyr, 
that dominates the optical spectrum, extended across 65 
arcsec {^^ 24 kpc), and covering the measurable extent of 
the galaxy. 

• YSP component. We have found a young stellar com- 
ponent, with ages tysp < 100 Myr, with varying contribu- 
tions across the galaxy, but particularly significant in the 
nuclear regions, where it is highly reddened. 

• Hidden Starburst. Our results reinforce the idea that 
there is prodigious ongoing star formation activity that 
is hidden at optical wavelengths, likely to be related to 
the radio supernov ae detected in the center of the galaxy 
l|Smith et al.l . Il998l ) , as well as the powerful mid to far-IR 
emission from this object. 



4 DISCUSSION 

4.1 Comparison with previous results 

IWilson et al. I l|2006t ) found evidence for a centrally located 
YSP (< 10 Myr) and an intermediate-age (~ 300 Myr) stel- 
lar population spread towards the north of the galaxy in 
their photometric study of the cluster populations in Arp 

* ISP of 700 Myr {E(B -V) = 0.2) + a YSP of 5 Myr {E(B - 
V) = 0.5) making a contribution of 80% and 20% of the total flux 
respectively 



220. These results are broadly consistent with the results 
presented in this paper. In addition, we find evidence for 
ISPs at all locations in the galaxy. 

There have been relatively few attempts in the past 
to study the stellar populati o ns in th e diff use hght of 
ULIRGs. ICanalizo fc StocktonI l|2000al lbl [20o3) found stel- 
lar ages ranging from currently ongoing starburst activ- 
ity to poststarburst ages of < 300 Myr in their sample of 
"transition QSOs" — objects that may represent a transi- 
tionary stage between ULIRGs and QSO. The ages found 
for the dominant stellar populations in Arp220 are gen- 
erally o ld er th an 300 Myr, although ICanalizo fc StocktonI 
(200Qa||3, 1 

200 ll ) used only one stellar component for their 

modelling. 

More recently, iRodrfguez Zaurm et al.l |2003) applied 
the same techniques used here to the ULIRG/radio galaxy 
PKS1345-I-12. Overall, the results found in their paper are 
consistent with the results found here for Arp 220, with a 
mixture of YSP and ISP components. However, precise age 
determination for the ISP component is more difficult in 
PKS1345-I-12 because of the presence of a strong OSP com- 
ponent. 

We are currently analyzing the optical spectra of a sam- 
ple of 40 ULIRGs (Rodriguez et al., 2007, in prep). We find 
that the modelling results are consistent with a combination 
of YSP -I- ISP components for the majority of the objects 
in our sample. However, we find that the YSPs found in 
these galaxies make a larger contribution than in the case 
of Arp220, and sometimes dominate the optical spectra. Fi- 
nally, similar results in terms of the mix of stellar popu- 
lations have recen tly been obtained i n det ail spectroscopic 
studies of LIRGs l|Gonzalez Delgadd . [2007h . 

4.2 The origin of the ISPs 

In order to use studies of the stellar populations to inves- 
tigate star forming histories in major galaxy merges it is 
important to understand the extent to which the stellar 
populations have been formed during the merger. A poten- 
tial issue is contamination by the young stellar populations 
present in the galaxy disks prior to the star of the merger. 
As a first approach to investigate the origin of the ISPs we 
used the CONFIT code and co mbined unreddened S a, Sb 
and Sc galaxy template spectra (jKinnev et al.l . Il996f ) with 
a stellar population with ages in the range 0.001 < age < 5 
Gyr and reddenings 0.0 < E{B -V) < 2.0. 

No good fits were found for combinations of Sa, Sb or 
Sc galaxy templates without a significant additional con- 
tribution of a YSP or an ISP component. However, the 
iKinnev et al.l l|l996l ) templates represent an average for sev- 
eral galaxies of the same morphological type, and each type 
encompasses a range of individual galaxy spectra, some of 
which deviate substantially from the averaged SED. There- 
fore, although the modelling results obtained using this com- 
bination suggest that the stellar populations detected in the 
optical for the ULIRG Arp220 have been formed during the 
merger event rather than captured, we cannot entirely rule 
out the idea that the ISPs detected in the galaxy represent 
the disrupted disks of one or more of the merging galaxies. 

Another way to address this issue is to compare the 
mass of YSPs in late type spiral galaxies, the likely progen- 
itor galaxies, with the values found here for the ISPs in Arp 
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Figure 4. 
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PA 160: Ap F: 0.5 Gyr ISP, E(B - V) = 0.4 + 0.005 Gyr, E(B - V) = 0.5 
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Figure 4. Detailed fits to the spectra extracted for some of the apertures sampling the different regions of the galaxy (see Figure 1 for 
the location of the extraction apertures). The green and blue lines are the ISP and YSP components respectively, while the red lines are 
the sum of both. It is clear from the figure that a variety of reddenings is required to fit the data for the different apertures. Also shown 
in the figure are detailed fits of the wavelength range 3700 A — 4300A (plots on left). The fluxes are presented in wavelength units. 
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STARLIGHT 
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0.2 
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0.2 
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0.2 
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14 
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Table 2. Modelhng results obtained using the STARLIGHT code with the Gonzalez Delgado et al. (2005) synthetic SEDs. Col (1): Slit 
PA. Col (2): Aperture label. Apertures D and E from PA 160 sample the nuclear regions of the galaxy Col (3): Width of the aperture 
(the centroid of Ap D is chosen to be the reference point for all the slit positions) . Col (4) : The intrinsic reddening value obtained for the 
best fitting models assuming the same extinction for all the stellar components. Cols (5), (6) and (7): The flux percentage at 4020A of 
each population, as defined in Section 3. 



220. The reason to focus only on YSPs in the parent galaxies 
is that, if the star formation in the disks is truncated during 
the merger, these are the populations which may evolve into 
ISPs detected today in Arp 220. Since our three slit positions 
do not cover the full extent of the object, we used the ap- 
parent magni t ude o f the galaxy in the V band, ray = 13.2 
l|Tavlor et al.1 . [200^, and applied the same techniques de- 
scribed in section 3.4 to account for possible slit losses. We 
obtain total masses in the range of 2.7x10"' -4.3xlO^°M^. 
for the ISPs in Arp220. To estimate the stellar mass content 
of lat e type spiral galaxie s, we modelled a Sc galaxy tem- 
plate ()Kinnev et al.l . ll996l ) combining a 12.5 Gyr OSP with 
a YSP. We find that we can model the template with an OSP 
plus YSPs of age tYSP< 20 Myr contributing 70 - 100% of 
the total flux in the normahsing bin (4420 A- 4500 A, cho- 
sen to be close to the B band wavelength range). We then 
assumed the same values for the contribution of each compo- 
nent to the median B-band monochro matic luminosities of 
Sc galaxies |Roberts fc Havnesl ll99J ). Since the templates 
are generated for 1 Mg , by scaling we were able to estimate 
the total mass and the mass contribution of each of the com- 
ponents. We estimated a range of 3.6x10^ -3.7x10^^Mq for 
the total stellar mass of the pop ulation of typical S c galax - 
ies in the large compilation of iRoberts fc Havnea l)l994h , 
of which only 1.2x10^ - 8.5 xlO^Mg is accounted by the 
YSP component. Considering only the top 25% of the most 



massive Sc galaxies in the lRoberts fc Havnesl (|l994h sample, 
the latter values increase by a factor of < 2. Therefore, it 
is possible that the ISPs detected in Arp220 represent the 
captured disks of one or more of the parent galaxies, but 
only if these galaxies are among the most massive late type 
spirals. 

In terms of gas content, IScoviUe et aP |l993) found a 
mass of ~ 3 x lO^^Mg for the H2 mass of Arp 220. How- 
ever, they used a standard Galactic CO-to-H2 ratio, which 
is likely to over estimate the molecula r gass mass M(H2) by a 
factor of three jSolomon et al.l . ll997h . On the other hand, an 
upper limit for the HI mas s value for Arp220 is 5.0xlO^M„ 
(|Mirabel fc Sandersl[l98^ ). Therefore, the total gas mass in 
the galaxy, M(H2 + HI), is of order of ~ 1.5 xlO^°M0. In 
comparison, the median value of the HI mas s cont ent of Sc 
galaxies is 0.8x10^°Mq (|Roberts fc Havnesl Il994l). Asum- 



ing a typical H2/HI mass ratio of 0.7 ( Young fc Kneze"^ 
Il989l ) for such galaxies, we find a gas mass content of M(H2 
-I- HI) ~ 1.4x10^°Mq for a typical Sc galaxy. Therefore, we 
conclude that it is possible for a merger of two typical spiral 
galaxies to account for the gas content estimated in Arp220 
provided that a large fraction of the ga s is not ejected dur- 
ing the merger (jPi Matteo et all 120071 ) or transformed into 
stars. However if a substantial proportion of the ISPs are 
formed in the merger, then the total amount of gas required 
is much larger, implying a merger between two Sc galaxies 
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CONFIT 



Width of Age of E{B - V) Age of E(B -V) %YSP 

the Apert ISP YSP of total 

(arcsec) (Gyr) (Gyr) flux 

(1) (2) (3) (4) (5) (6) (7) (8) 



PA 160 Ap A -12.6 to -8.2 0.6 - 0.7 0.0 - 0.2 < 0.1 0.0-0.2 <5 

Ap B -8.2 to -5.4 0.5 - 0.7 0.0 - 0.2 < 0.1 0.0-0.3 < 20 

Ap C -5.4 to -1.8 0.5 - 0.7 0.3 - 0.5 < 0.1 0.0-1.0 5-15 

Ap D -1.8 to 1.8 0.5 - 0.9 0.0 - 1.0 < 0.08 0.9 - 1.8 30 - 63 

Ap E 1.8 to 4.6 0.5 - 0.9 0.2 - 0.6 < 0.06 0.9 - 1.6 22 - 47 

Ap F 4.6 to 8.6 0.5 - 0.7 0.2 - 0.6 < 0.08 0.0 - 1.5 5 - 36 

Ap G 8.6 to 22.6 0.5 - 0.7 0.0 - 0.3 < 0.07 0.0 - 0.7 5 - 20 

APtotal -24.4 to 28.4 0.5 - 0.7 0.2 - 0.6 < 0.1 0.0-1.5 5-35 



PA 75 Ap A -26.6 to -8.6 0.5 - 0.7 0.2 - 0.5 < 0.1 0.0 - 2.0 < 20 

Ap B -8.6 to -4.2 0.5 - 0.7 0.2 - 0.6 < 0.1 0.0-1.0 5-25 

Ap C -4.2 to -1.4 0.5 - 0.7 0.2 - 0.6 < 0.05 0.0 - 1.2 5 - 30 

Ap D -1.4 to 1.4 0.5 - 0.7 0.2 - 0.6 < 0.06 0.0 - 1.5 5 - 25 

Ap E 1.4 to 5.0 0.5 - 0.7 0.2 - 0.5 < 0.05 0.0 - 2.0 < 30 

Ap F 5.0 to 10.2 0.6 - 0.9 0.1 - 0.4 < 0.1 0.0-1.2 < 25 

Ap G 10.2 to 34.6 0.7 - 0.9 0.0 - 0.3 < 0.1 0.0-0.8 < 25 

APtotal -29.4 to 35.8 0.5 - 0.7 0.2 - 0.5 < 0.1 0.0-1.2 5-25 



PA75* Ap A -19.2 to -12.0 0.5 - 0.7 0.0 - 0.3 < 0.06 0.0 - 0.6 10 - 35 

Ap B -12.0 to 5.2 0.5 - 0.7 0.0 - 0.3 < 0.1 0.0-1.5 5-25 

Ap C -5.2 to -2.8 0.5 - 0.7 0.0 - 0.3 < 0.1 0.0-1.5 < 40 

Ap D -2.8 to 2.8 0.5 - 0.8 0.0 - 0.3 0.06 0.0 - 1.5 < 40 

Ap E 2.8 to 17.6 0.6 - 0.9 0.0 - 0.3 < 0.06 0.0 - 1.0 < 25 

Ap F 17.6 to 38.4 0.6 - 0.9 0.0 - 0.3 < 0.1 0.0-2.0 < 15 

APtotal -23.2 to 45.6 0.5 - 0.8 0.0 - 0.3 < 0.06 0.0 - 1.5 5 - 30 



Table 3. Modelling results obtained using the CONFIT code with the Bruzual & Chariot (2003) synthetic SEDs. Col (1): Slit PA. Col 
(2): Aperture label. Apertures D and E from PA 160 sample the nuclear regions of the galaxy Col (3): Width of the aperture (the centroid 
of Ap D is chosen to be the reference point for all the slit positions). Col (4); Range of ages for the ISP component. Col (5): Range of 
intrinsic E{B — V) values for the ISP component. Col (6): Upper limits for the ages of the YSP component. Col (7): Same as column 
5 for the YSP component, col (8): Upper limits for the contribution in flux of the YSP to the model in the normalising bin (4700A — 
4800A), for those cases for which an ISP can model the data with a negligible contribution of a young component. 



in th e upper 25% of the HI mass range (|Roberts fc Havnesl 
1 19941 ). 



4.3 Comparison with models 

Arp220 is a double nucleus system with a separation of 0.95 
arcsec, corresponding to ~345 pc, that exhibits tidal tails in 
the optical. The morphology of the galaxy suggests that it 
is in the final stages of the merger event, just before the two 
nuclei coalesce. In this section we will compare the results 
found in this paper with the prediction of the models for the 
star formation activity at such a stage of a merger event. 

In ge neral, simulations pr e dict t w o epochs of starburst 
activ i ty (iMihos fc Hernguistl , Il996l : iBarnes fc Hernguistl 
1 19961 : ISpringel et all, booj ) in major gas-rich mergers: the 
first ocurring just after the first encounter, and the sec- 
ond, more intense, episode towards the end of the merger, 
when the nuclei coalesce. However, both the time lag and 
the relative intensity of the peaks of starburst activity dur- 
ing the merger event depend on several factors: the pres- 
ence of bulges, feedback effects, gas content and orbital 
geometry. For example, the presence of a bulge acts as 
a stabilizer of the gas against infiows and formation of 
bar structures, allowing stronger starburst activity t o wards 
the end of the merger event (|Mihos fc Hernguistl 1 19961 : 



iBarnes fc Hernguistll 19961 ). On the other hand, AGN feed- 
back effects (e.g. quasar-driven winds) disrupt the gas sur- 
rounding th^_black_Jiolej_a£ting against the star formation 
activity (|Springel et all . boOSi ) . 

In the case of Arp220, we derive star formation rates of 
~ 260 Mq yr~^, consistent with a starburst being respon- 
sible for the IR luminosity of the galaxy and there is no 
evidence for AGN activity (jGenzel et al.l. 19981). B ased on 
the results of the simulations I Springel et al.l . 2005h . highly 
gas-rich disc galaxies must be involved in the merger event 
to account for such high star formation rates. In this con- 
text, the results found in the previous section are consistent 
with the model predictions. It is possible that the dominant 
ISPs detected in Arp220 are associated with the first en- 
hancement in star formation activity that occurs in the early 
stages of the merger, coinciding with the first encounter be- 
tween the merging nuclei. However, because of the potential 
contamination from stars in the disrupted disks of the merg- 
ing galaxies (section 4.2), the fraction of ISPs formed in the 
merger is not known with any accuracy. On the other hand, 
the YSPs detected in the center of the galaxy, as well as 
the ULIRG activity, are likely to be related to the major 
enhancement of the activity as the two nuclei coalesce in 
the final stages of the merger. It is notable that the ages 
of the YSP (tysp < 0.1 Gyr) are consistent with the 0.1 
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Gyr timescale of the final starburst predicted by the merger 
simulations. 



5 CONCLUSIONS 

The results presented in this paper demonstrate the utility 
of spectroscopic studies of the diffuse light for investigating 
the evolution of the stellar populations in ULIRGs. 

We have found evidence for a complex star forming 
history in the ULIRG Arp 220, with at least two distinct 
episodes of starburst activity. The results are consistent with 
an intermediate-age stellar population of age 0.1 < tjsp < 
0.9 Gyr, dominating the optical emission throughout the 
body of the galaxy. We have also detected young (< 100 
Myr) stellar populations (YSPs) located in the central ~ 
2.5 kpc that contribute significantly to the optical emission. 
This latter population is likely to be related to the last en- 
hancement of the activity towards the end of the merger 
event. 

On the other hand, we have estimated that the bolo- 
metric luminosity accounted by the optically visible stel- 
lar populations and compared it with the IR-luminosity of 
the galaxy. We conclude that we are not directly detect- 
ing the bulk of the starburst activity responsible for the 
ULIRG phenomenon in Arp220, whic h is likely to bere- 
lated to the radio supernovae found bv lSmith et aD (|l998h . 
representing another "extra" episode of star formation ac- 
tivity apart from the two suggested by our modelling re- 
sults. It is clear that our results for Arp220 are consistent 
with the emerging evidence for a complex and multimodal 
star formation a c tivity in merging s ystems ( The An t ennae : 
IWhitmore eraP (Il999h. NGC7252: iMaraston et aD (|200lh . 
PKS1345-H2: lRodn'guez Zaurm et al.1 l|2007l )). 

In terms of reddening, while the presence of highly red- 
dened stellar component is unlikely in the extended regions 
of the galaxy, high reddening values are required to match 
the data for the apertures sampling the nuclear regions, co- 
inciding with dust lanes in the HST images. Our results show 
the importance of accounting for reddening when trying to 
model stellar populations in star forming galaxies. 

Finally, we note that the modelling results for the 
three large apertures, representing the integrated light of 
the galaxy, are consistent with those of the detailed study 
based on smaller apertures. This gives confidence to studies 
of ULIRGs at medium and high redshifts, suggesting that 
age measurements based on relatively large aperture studies 
of the diffuse light at such redshifts, are representative of 
the stellar populations of the galaxies as a whole. 
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